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Abstract — A new mechanism describing the formation of protogalaxies is proposed, which is based on the sec- 
ond-order phase transition in the inflation stage and the domain wall formation upon the end of inflation. The 
presence of closed domain walls with the size markedly exceeding the cosmological horizon at the instant of 
their formation and the wall collapse in the postinflation epoch (when the wall size becomes comparable with 
the cosmological horizon), which leads to the formation of massive black hole clusters that can serve as nuclei 
for the future galaxies. The black hole mass distributions obtained do not contradict the available experimental 
data. The number of black holes with M ~ 100 solar masses (Mq) and above is comparable with the number of 
Galaxies in the visible Universe. Development of the proposed approach gives ground for a principally new sce- 
nario of the galaxy formation in the model of hot Universe. © 2001 MAIK "Nauka/Interperiodica". 



1. INTRODUCTION 

In the past decade, investigations into the nature of 
active galactic nuclei exhibited a considerable progress. 
Now there is virtually no doubt that the centers of gal- 
axies contain massive black holes [1]. It is the existence 
of black holes with masses on the order of 10^-10** Mq 
in the galactic nuclei and the accretion of matter onto 
these holes that is believed to account for the physical 
nature of their activity. A possible explanation of the 
formation of such supermassive black holes assumes 
the collapse of a large number of stars caused by their 
high concentration at the galaxy center. However, the 
mechanism of the galactic nuclei formation is still 
unclear. According to Veilleux [2], there are serious 
grounds to believe that the formation of stars and gal- 
axies proceeded simultaneously. Stiavelli [3] consid- 
ered a model of the galaxy formation around a massive 
black hole and presented arguments in favor of this 
model (see also [4]). Each of the two approaches has 
certain advantages, while being not free of drawbacks. 

The problem of the possible "primordial" black hole 
(PBH) formation is still open. In contrast to the case of 
"secondary" black holes, which are related to the evo- 
lution of stars and stellar systems, there is no convinc- 
ing astronomic evidence for the existence of PBHs. 
Restrictions posed by the astronomic data on the PBH 
concentration offer a unique source of information on 
the inhomogeneity of early Universe [5] and on the 
physical processes accounting for this inhomogeneity 



[6]. Generally speaking, the PBH mass may be arbi- 
trary, ranging from the Planck value (or even below [7]) 
up to a mass contained within the contemporary cos- 
mological horizon. However, in most cases the astro- 
physical effects related to the presence of PBH are 
restricted to masses much lower than the solar mass. 
The reason is that the mechanism of PBH formation is 
usually related to the development of inhomogeneities 
bounded by the cosmological horizon. The data of 
observations concerning the distribution of light ele- 
ments and the spectrum of cosmic microwave radiation 
pose very rigid restrictions on the magnitude of inho- 
mogeneities existing in the pregalactic stage following 
the first second of expansion of the Universe. Thus, 
realistic mechanisms of the PBH formation have to be 
apparently related to very early (f <l 1 s) stages of evo- 
lution of the Universe — when the mass contained 
within the cosmological horizon and limiting the possi- 
ble PBH mass was significantly lower than the solar 
mass. Nevertheless, the actively discussed possibility 
of a genetic relationship between quasars and active 
galactic nuclei on the one hand and the existence of 
PBHs with much greater masses on the other hand [8] 
becomes a subject for detailed investigations [9-1 1]. 

Below we will consider a new mechanism describ- 
ing the early formation of PBHs, which serve as the 
nucleation centers in the subsequent formation of gal- 
axies. This mechanism may prove to be free from dis- 
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advantages inherent in the models based on the concept 
of a single PBH being a nucleus of the future galaxy. 

Previously [11 J we proposed a new mechanism of 
the PBH formation that opens possibihty of the massive 
black hole formation in the early Universe. The mech- 
anism is based on the possibility that black holes are 
formed as a result of a collapse of closed walls formed 
during a second-order phase transition. The masses of 
such black holes may vary within broad hmits, up to a 
level on the order of 10^ Mq. 

Let us assume that a potential of the field system 
possesses at least two different vacuum states. Then 
there are two possible distributions of these states in the 
early Universe. The first possibility is that the Universe 
contains approximately equal numbers of both states, 
which is typical of a temperature- controlled phase 
transition under usual conditions. The alternative possi- 
bility corresponds to the case when the two vacuum 
states form with different probabilities. In this case, 
there appear islands of less probable vacuum state sur- 
rounded by the sea of another, more probable vacuum 
state. As was recently demonstrated [12], an important 
condition for this distribution is the existence of valleys 
in the scalar field potential during inflation. Then the 
background de Sitter fluctuations in this massless scalar 
field lead to the formation of islands representing one 
vacuum in the sea of another vacuum, this phase transi- 
tion takes place after the end of inflation in the Fried- 
mann-Robertson-Walker Universe. After the phase 
transition, the two vacuum states are separated by a 
wall; the size of this wall may be significantly greater 
as compared to the cosmological horizon at that period 
of time. At some instant after crossing the horizon, the 
walls begin to contract because of the surface tension. 
As a result, provided that friction is absent and the wall 
does not radiate a considerable part of its energy in the 
form of scalar waves, almost all energy of this closed 
wall may be concentrated within a small volume inside 
the gravitational radius. This is a necessary condition of 
the black hole formation. 

The mass spectrum of black holes formed by this 
mechanism depends on parameters of the scalar field 
potential determining the direction and size of the 
potential valley during inflation and the postinflation 
phase transition. Although we deal here with the so- 
called pseudo-Nambu-Goldstone field, the proposed 
mechanism possesses a sufficiently general character. 

The presence of massive PBHs is a new factor in the 
development of gravitational instability in the sur- 
rounding matter and may serve a base for new scenarios 
of the formation and evolution of galaxies. 

2. THE FORMATION OF CLOSED WALLS 
IN A COMPLEX FIELD 

Now we will describe a mechanism accounting for 
the appearance of massive walls with the size markedly 



greater than the horizon at the end of inflation. Let us 
consider a complex scalar field with the potential 

y((p) = uW\^-fi2y, (1) 

where (p = re*. We assume the mass of the radial field 
component r to be sufficiently large, so that the com- 
plex field would occur in the ground state even before 
the end of inflation. Since the minimum of potential (1) 
is degenerate, the field has the form 

(p-(//72y', 

with the phase acquiring the meaning of a massless 
field. 

For the following considerations, it should be noted 
that, using expression (1) in the inflation period, we 
ignored the term 

8y(e) = A''(i-cose). (2) 

reflecting the contribution of instanton effects to the 
Lagrangian renormaUzation. Since the parameter A 
appears as a result of the normahzation, its value cannot 
be large and we may quite reasonably assume that A < 
H, f. The omitted term (2) begins to play a significant 
role in the postinflation stage, when the Hubble con- 
stant sharply decreases with time {H = l/2t during the 
radiation dominated epoch). 

Let us assume that a part of the Universe occurring 
inside the contemporary horizon was formed N^ 
e-folds before the end of inflation. As was demon- 
strated in [13], the quantum field fluctuations during 
inflation are rapidly transformed into a classical field 
component, while the massless field values in the 
neighboring causahty-disconnected intervals with the 
size H'^ differ on the average by 

56 = H/2nf (3) 

after a single e-fold. In the next time step At - (i.e., 
during the next e-fold) each causahty-connected 
domain is divided into ~e^ causality-disconnected sub- 
domains; the phase in each of the new domains differs 
by -89 from that in the preceding step. Thus, more and 
more domains appear with time in which the phase dif- 
fers significantly from the initial value. More precisely, 
the probability of finding the phase is [14-16] 

P,e,«, = ^exp{-t^}. 

where N is the number of e-folds remaining to the end 
of the infiation period, 9;; is the random phase value at 
the instant of formation of the causality-connected 
domain corresponding to a visible part of the contem- 
porary Universe. Without a loss of generality, we may 
select < Gf; < :i. Below we wiU demonstrate that a par- 
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ticular value of the initial phase significantly affects 
evolution of the Universe in the postinflation epoch. 

The probability of finding a certain phase obeys the 
Gaussian distribution (4) and, hence, the phase aver- 
aged over the entire space equals the random initial 
phase 9y. A principally important point is the appear- 
ance of domains with the phases 6 > 7i. Appearing only 
after a certain period of time during which the Universe 
exhibited exponential expansion, these domains turn 
out to be surrounded by a space with the phase < :i. 
These very domains lead in the following to the forma- 
tion of large-scale structures. Note that the phase fluc- 
tuations during the first e-folds may, generally speak- 
ing, transform eventually into fluctuations of the cos- 
mic microwave radiation, which will lead to imposing 
restrictions on the scaling parameter/. This difficulty 
can be avoided by taking into account interaction of the 
field (p with the inflaton field (i.e., by making parameter 
/a variable). 

Initially, the potential (1) possessed a U{1) symme- 
try and the phase corresponded to a massless scalar 
field. Owing to the term (2), the symmetry is broken 
after the end of the inflation period: the potential of 
the 9 field acquires minima at the points = 0, +2ji, 
±4n, and the field acquires the mass = 2//A^. 
According to the classical equation of motion. 



+ 3//0 + 



d5V 



= 0, 



(5) 



the phase performs decaying oscillations about the 
potential minimum, the initial values being different in 
various space domains. Moreover, domains with the 
initial phase Jt < 9 < 27t perform oscillations about the 
potential minimum at 9„i„ = 2jt, whereas the phase in 
the surrounding space tends to a minimum at the point 
9min = 0. Upon ceasing of the decaying phase oscilla- 
tions, the system contains domains characterized by the 
phase 0^n = 2n surrounded by the space with 0„j„ = 0. 
Apparently, on moving in any direction from inside to 
outside of the domain, we will unavoidably pass trough 
a point where 9 = Jt because the phase varies continu- 
ously. This imphes that a closed surface must exist 
which is characterized by the phase Q^^i = n. The size 
of this surface depends on the moment of domain for- 
mation in the inflation period, while the shape of the 
surface may be arbitrary. A principal point for the sub- 
sequent considerations is that the surface is closed. 

Thus, we obtained a field configuration connecting 
various vacuum states of the potential (2). A rigorous 
classical solution of this problem possessing a transla- 
tional symmetry in the two directions in space (flat 
wall) is well known [17]: 



%Ax-Xo) = -4arctan 



exp^— 



(6) 



where d is the wall thickness. Since the thickness of a 
closed wall is related to microscopic parameters of the 



theory, whereas the characteristic wall size is a priori 
unlimited, expression (6) is applicable to within a suf- 
ficiently high accuracy in the case under consideration. 
As can be readily shown, the wall possesses an energy 
concentrated where the phase is 9 = tt [17J. Thus, we 
obtained a mechanism providing the formation of 
domains surrounded by closed walls. The surface energy 
of a wall depends on the Lagrangian parameters, while 
the wall size id determined by the time of crossing the 
phase value equal to n during the inflation period. 



3. DOMAIN WALL DECELERATION DURING 
MOTION THROUGH A PLASMA 

The first- and second-order phase transitions lead to 
the formation of a field walls separating one vacuum of 
this field from another One of these mechanisms was 
described in the preceding section. In turn, the walls are 
moving at a subluminal velocity and interact with the 
surrounding plasma. Depending on the character of this 
interaction and the shape of the field potential, there are 
two possible situations. In the first case, the plasma par- 
ticles pass through the waU, faUing into a different vac- 
uum and acquiring a certain mass. This situation corre- 
sponds to an electroweak phase interaction [18], 
whereby the corresponding Higgs field is responsible 
for a mechanism of the fermion mass production. In the 
opposite case, the particle mass does not change upon 
going from one to another vacuum (an example is 
offered by the case of interaction with an axion wall). 
In the former case, the interaction with the medium 
leads to a significant retardation of the domain wall, 
while in the latter case, the wall are virtually transpar- 
ent for the medium provided that the parameters are 
given reasonable values. 

All considerations are convenientiy conducted in 
the resting wall frame. The probability of a particle 
scattering from the plane resting wall is 



dw = dn{k)2n8{t-t')\M\- 



2ey(23i) e' 



(7) 



where dn(k) is the distribution of incident particles with 
respect to momentum, M is the matrix element for the 
particle transition from a state with the energy e and 
momentum k to the state with the energy e' and momen- 
tum k' upon interaction with the potential U = U{z) 
describing the plane wall. The pressure produced by 
incident particles upon the wall is related to the rate of 
their momentum transfer to the wall 



(8) 



where S is the wall area. Let us select the Lagrangian of 
the particle-wall interaction in the foUowing form rep- 
resenting a classical configuration of the complex field 
phase: 

Lint = Kdfi{z)J„ = xi/Y^¥. K = i/f. (9) 
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Calculating a matrix element for the particle scatter- 
ing from the waU with the transition from the initial 
momentum k to final momentum k' 



M = {k'\^L;^,(fx\k) 



we obtain 



|M|' = S{4nfKS8^^\q«)k^ 



z 2 * 

cosh {k dn) 



(10) 



(11) 



In deriving formula (11), we took into account that 
the laws of the energy-momentum conservation lead to 
the following relationships 

k[ = ±k^; = /c|| - /c|| = 0, 

according to which a nonzero contribution to the pres- 
sure is only due to the reflected particles with k[ = -k^. 
Therefore, the pressure of incident particles upon the 
wall can be written as 



4 2f k, 

= — 

n ■' cosh (nkJ) 



{k,-k'ME-E') 



(12) 



^Hk,-kl/-fl'^. 



Let us determine a distribution of the incident parti- 
cles with respect to the transverse momentum dn{k). In 
the resting plasma frame. 



dn^ikg) = Cexp 



E„{ko)]d\oV 



1 



(13) 



Here and below, the subscript denotes quantities 
determined in the resting plasma frame. Assuming the 
plasma temperature T to be significantly greater as 
compared to the fermion masses and normaUzing it to 
the total particle density 

«,ot = iVte*)r', iV(^*)-5, 

we obtain J = 5/871. In addition, it is evident that 

dn{k) = dn„{kf^), (14) 

where the incident particle momentum in the resting 
wall and plasma frames (in the latter frame, the waU 
moves at a velocity v) are related as 

^0|| = ^11' 

(15) 

Eo = yivk, + e), 
1 



Integrating the pressure (12) with resect to momen- 
tum of the incident particle, we obtain 



,32Ck^ cd^k, . , k^ 



p = sgn(A:,)— -jY — (E+vA:,)- 

(271)' J 8 cosh (Tik,d) 



(16) 



X exp-^ -y 



e+ vk. 



This formula was derived with an allowance of the 
Lorentz invariance of the phase volume d^k'lt'. Numer- 
ical calculation of the integral in (16) presents no diffi- 
culties, but we are interested in analytically estimating 
the pressure produced by the medium upon the wall. 
For this purpose, note that the walls are formed at tem- 
peratures r ~ A and the wall thickness is = //2A^. 
Therefore, there is a large parameter 

™-//A > 1, 

using which we may obtain a sufficiently reliable esti- 
mate of the integral. According to (16), the most effec- 
tive scattering takes place for an incident particle 
momentum of 

k,~\lTid<T. 
At the same time, it is evident that 

k^^~t~rr>k,. 

Using these relationships, we may estimate the inte- 
gral in (16). A final expression for the pressure pro- 
duced by the surrounding medium upon the domain 
wall is as follows: 



20r y 

7t^ 



(17) 



4. CONDITIONS FOR THE PRIMORDIAL 
BLACK HOLE FORMATION 

After heating of the Universe, the evolution of 
domains formed with the phase > 7l and sharply 
increased in volume during the inflation period pro- 
ceeds on the background of the Friedmann expansion 
and is described by the relativistic equation of state. 
First, an equilibrium state with the "vacuum" phase = 
271 inside the domain and the = phase outside is 
established at A. A closed wall corresponding to the 
phase 9 = 71 is formed in the transition region with a 
width of ~\lm ~f/A^, which separates the domain from 
the surrounding space. The surface energy density on 
the wafl amounts to ~//A^. 

It must be noted that the process of establishing of 
the equihbrium ("vacuum") phase values may acquire 
a protracted character. If the stage of coherent phase 
osciflations about the equihbrium values is sufficiently 
long, the energy density of these oscillations may 
become dominating and determine the dust period of 
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expansion. Let us consider the factors influencing the 
cosmological evolution of such a wall. 

1. First, note that the domain size immediately after 
the end of inflation markedly exceeds the horizon size 
in the Friedmann expansion stage. The overall contrac- 
tion of the closed wall may begin only when the hori- 
zon size R), will be equal to the domain size. Up to this 
moment, the characteristic domain size increases with 
the expanding Universe because we assumed that the 
existing field contribution to the total energy-momen- 
tum tensor is small as compared to that of the inflaton 
field. Accordingly, the field gives also a small contribu- 
tion to the total energy density of the Universe upon 
heating, when the energy density of infiaton transforms 
into the energy density of relativistic particles. Evi- 
dently, internal stresses developed in the wall after 
crossing the horizon initiate processes tending to mini- 
mize the wall surface. This implies that the wall tends, 
first, to acquire a spherical shape and, second, to con- 
tract toward the center. For simplicity, below we will 
consider the motion of closed spherical walls. 

2. Since the energy of the surrounding plasma rap- 
idly decreases, the wall energy may become at a certain 
time instant comparable with the energy of the sur- 
rounding medium. Simultaneously, the domain sepa- 
rates from the general expansion process and its racfius 

may become smaller than Rf,. 

3. The wall energy is proportional to its area at the 

instant of crossing the horizon. By the moment of max- 
imum contraction, this energy is virtually completely 
converted into the kinetic energy. Should the wall by 
this moment be locahzed within the hmits of the gravi- 
tational radius, a PBH is formed. 

4. Contracting under the action of internal forces, 
the wall moves through the surrounding plasma. The 
resulting force of friction may, under certain condi- 
tions, become significant and lead to a uniform (nonac- 
celerated) contraction of the wall. In this case, the 
potential energy of the wall is dissipated in the sur- 
rounding medium. Only when the wall would decrease 
to a certain small size Rf, the internal forces propor- 
tional to the surface curvature wiU dominate and the 
wall will again contract with acceleration to supply a 
necessary energy to the center, sufficient to form a 
PBH. 

The above considerations show that the energy con- 
centrated in the course of waU contraction can be deter- 
mined using the condition 



£-4jt7?o, R = mm{Ri„R^^,Rf). 



(18) 



where a = AIs?f is the surface energy density of the 
wall. A condition of the PBH formation is that 



^min ~d<rg = lE/rripi. 



(19) 



It is assumed that the spherical waU contracts until 
reaching a size on the order of the wall thickness. 



Let us determine the values of R,,, R„ and Rf for a 
system with Lagrangian (1). Consider a domain with a 
certain phase appearing e-folds before the end of 
infiation. The domain size at the end of infiation period is 



l=H e , 



(20) 



where //is the Hubble constant at that time instant. By 
the moment of crossing the horizon, the domain will 
expand to acquire the characteristic size 



2HN' 



(21) 



Below we assume that a visible part of the Universe is 
formed N=Nij = 60 e-folds before the end of infiation. 

The second characteristic size R„ is determined 

from the condition that the wall energy (18) is equal to 
the energy of plasma contained in the domain bounded 
by the closed wall: 

Ey - 

where 

^ 30^ 

is the plasma density during the radiation dominated 
epoch. Taking into account that a = 4AV/, we obtain the 
critical wall size corresponding to the domain separat- 
ing from the general expansion: 



^ = ^crit = — 
r 



_ 3a _ 360aV _ aV 



(22) 



As is known, the temperature in the Robertson- Walker 
Universe during the radiation dominated epoch varies 
with time as 



T = 



45 



32ng'> 



(23) 



Taking into account that an increase in the wall radius 
up to the instant of separation from the general expan- 
sion is proportional to the scaling factor 

Rit) = hjtit,, (24) 
we arrive at the desired relationship for R„{N): 

R„ - i2'ncr[^) ^-exp(f ]. (25) 

An expression for the characteristic wall size (radius)/?^ 
above which the friction is significant can be obtained 
by equating the pressure developed by the internal 
forces Pint = 2a/Rfto that produced by the surrounding 
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medium on the moving wall. Using relationship (17) 
for the latter pressure, we obtain 



% cd 
10k J 



(26) 



The above considerations do not take into account the 
effect of a gravity field on the wall dynamics. There- 
fore, the obtained relationships are valid provided that 
the initial wall size is much greater than the gravita- 
tional radius. Generally speaking, the gravitational 
radius may be comparable with (or even exceed) the 
wall size for a sufficiently large domain size. In this 
study, we deal with smaller domains for which the 
intrinsic gravity field does not affect the wall evolution. 
Now we proceed to the study of PBH cluster formation 
in the early evolution stage of the Universe. 

5. CORRELATIONS 
IN THE BLACK HOLE DISTRIBUTION 

Previously [11] we have studied a new process 
involving the formation of uncorrected PBHs in the 
Universe. It was demonstrated that a model with rea- 
sonably selected parameters readily provides for the 
formation of 10" massive (10^°- 10"° g each) black 
holes, which is equal precisely equal to the number of 
Galaxies in the visible Universe. In that analysis, we 
did not took into account correlations (inherent in this 
mechanism) between the formation of a massive black 
hole and the appearance of smaller back holes sur- 
rounding the former one. This correlation is related pri- 
marily to certain features of the above-discussed pro- 
cess of the formation of domains with the phases > ;i. 
Apparently, the appearance of such domains creates 
prerequisites for the formation of new smaller domains 
inside. 

Let us estimate the mass distribution of these daugh- 
ter domains. Consider a region with a size on the order 
of H'^ and a phase within %<%<% + 8 (where 6 = 
H/2nfis the average phase jump during the time 
period) formed during the inflation period as a result of 
fluctuation in a certain region of space with the phase 
< 31. During the next e-fold, this space domain will 
separate into subdomains and some of these will 
acquire a phase 0i in the interval 31 - 5 < 0i < 7t. Upon 
the subsequent phase transition, these domains will be 
separated by walls from the external region. Similar 
transitions, with crossing the phase = Ji in the reverse 
direction will take place in each subdomain during the 
next e-fold. Thus, a structure of the fractal type appears 
which reproduces itself in each time step on a decreas- 
ing scale. 

Let ^ denote the number of subdomains formed in 
each step, around which a wall may form with time. 
Apparently, this value obeys the inequality 1 < ^ e^. 
In tiie subsequent estimates, we will assume that C, ~ 2-3. 
Since each causality-connected domain touches 



approximately six neighboring domains, we can hardly 
expect C, to be greater for a total number of = 20. 
The mass of the future black hole (if this would actually 
form) is determined by the area of a closed surface with 
the phase = The ratio of areas of the initial (mother) 
and daughter domains is readily estimated: the initial 
area after a single e- fold is 

So = e^H-\ 
and the daughter subdomain area is 
Si - H-\ 

Therefore, the ratio of masses of the black holes 
belonging to two sequential generations is 

M/Mj,,^S/Sf,,^e\ (27) 

for their relative number assumed to be 

Nj,i/N^ = t;. (28) 

As is readily seen, the number and mass of black holes 
appearing upon the y'th e-fold after the initial domain 
formation are determined by parameters of the largest 
black hole genetically related to the primary domain in 
which the phase originally exceeded n. It is evident that 



(29) 



Excluding j from these relationships, we obtain the 
desired black hole mass distribution in a cluster: 



NJM)^{Mo/M) 



(l/2)ln5 



(30) 



The total mass of the cluster can be expressed through 
the mass Mq of the largest initial black hole: 



M,„, = Mo + CMi + C'M2+... = Mo 
+ Ce-X + (Ce"VMo+... = Mo[l-C/eY'. 



(31) 



As is seen, the total mass of the black hole cluster is 
only one and a half to two times greater than the largest 
initial black hole mass. The number of daughter black 
holes depends on the factors considered in the next sec- 
tion. 



6. DISCUSSION OF RESULTS 

In the preceding sections, we considered only the 
principal possibility of the formation of domain walls 
connecting adjacent vacuum states. We have used the 
formulas derived above to estimate efficiency of the 
proposed mechanism of the black hole cluster forma- 
tion. The numerical calculations were performed for 
the following values of parameters (which are consis- 
tent with the observed anisotropy in the cosmic micro- 
wave radiation): the Hubble constant at the end of infla- 
tion, H = 10'^ GeV; Lagrangian parameters,/= 1.77// 
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Fig. 1. Plots of the characteristic domain size (incUned soUd 
Une) and gravitational radius (dashed Une) versus the 
domain mass. Horizontal Une indicates the wall thickness. 
See the text for numerical values of parameters. 

and A = 5 Gev. The initial phase, at which the visible 
part of the Universe is formed by the time t^ ~ 60//"' to 
the end of inflation, controls the number of domains 
and, accordingly, the number of closed walls formed in 
the postinflation stage. This random value, not related 
to the Lagrangian parameters, must be selected taking 
into account the data of observations on the abundance 
of black holes in the Universe. We will use the numeri- 
cal value 6y = O.OSJt, which ensures a sufficiently large 
number of massive black holes, while the presence of 
numerous smaller black holes does not contradict 
experimental restrictions. 

Figure 1 shows the results of numerical calculations 
constructed in the logarithmic scale. The bottom hori- 
zontal line shows the wall thickness. As is seen, the 
condition of wall existence (the characteristic domain 
size must exceed the wall thickness) is fulfilled for the 
domains with masses exceeding 10'^ g. The domains of 
lower energies possess (for the parameters selected) 
dimensions below the wall thickness. This implies that 
the wall formation is impossible and the domain exhib- 
its only fluctuations in the energy density. During con- 
traction, the wall energy is eventually completely con- 
verted into radiation. 

The necessary condition for a black hole formation 
as a result of the domain wall collapse is that the grav- 
itational radius of the domain must be greater than the 
wall thickness. As is seen from fig. 1, this condition is 
fulfilled for the black holes with masses >102^ g. There- 
fore, the proposed mechanism of the black hole forma- 
tion leads to a nontrivial situation: massive PBHs exist 
at a compete absence of the black holes possessing 
masses below this threshold. Note that most significant 
observational restrictions concerning the PBH abun- 
dance refer to the mass region -10'' g (which cannot 
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Fig. 2. Primordial black hole mass distribution in the Uni- 
verse. 



form for the parameters selected). Since the gravita- 
tional radius is proportional to the wall surface area, the 
plots corresponding to the domain size and its gravita- 
tional radius must intersect. This intersection actually 
takes pace at a wall mass of 10^' g. For greater masses, 
the gravity effects have to be taken from the very begin- 
ning, which will limit from above the possible PBH 
masses. According to formula (26), the friction 
becomes significant only for supermassive walls, the 
number of which is negligibly small. 

Figure 2 shows the PBH mass distribution calcu- 
lated for the selected parameters (see also the discus- 
sion in [1 1]). As is seen, the PBH masses faU within the 
range from 102^ to 10^' g. The initial phase Q^ was 
selected so as to provide that the number of massive 
PBHs (-10^^ g) was equal to the number of Galaxies in 
the visible part of the Universe. The total mass of black 
holes amounts to -1% of the contemporary baryonic 
contribution. 

The results of calculations are sensitive to changes 
in the parameter A and the initial phase Q^. As the A 
value decreases to ~l GeV, still greater PBHs appear 
with a mass of up to -10'"' g. A change in the initial 
phase leads to sharp variations in the total number of 
black holes. As was shown in Section 5, each domain 
generates a family of subdomains in the close vicinity. 
A total mass of such a cluster is only 1.5-2 times that 
of the largest initial black hole in this space region. 
Thus, our calculations confirm the possibility of forma- 
tion of the clusters of massive PBHs (-100 M© and 
above) in the earhest stages of evolution of the Uni- 
verse at a temperature of 1-10 GeV. These clusters rep- 
resent stable energy density fluctuations around which 
increased baryonic density may concentrate in the sub- 
sequent stages, followed by the evolution into Galaxies. 
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1. CONCLUSION 

This paper proposes a new mechanism for the for- 
mation of protogalaxies, which is based on the cosmo- 
logical inferences of the elementary particle models 
predicting nonequilibriumsecond- order phase transi- 
tion in the inflation stage period and the domain wall 
formation upon the end of inflation. The presence of 
closed domain walls with the size markedly exceeding 
the cosmological horizon at the instant of their forma- 
tion leads to the wall collapse in the postinflation epoch 
(when the wall size becomes comparable with the cos- 
mological horizon), which results in the formation of 
massive black hole clusters that can serve as nuclei for 
the future galaxies. The black hole mass distributions 
are calculated which do not contradict the available 
experimental data. The number of black holes with M ~ 
100 Mq and above is comparable with the number of 
Galaxies in the visible Universe. A mechanism of 
deceleration of the wall motion is considered and it is 
shown that this process may affect only the dynamics of 
collapse of supermassive walls. 

Development of the proposed approach gives 
ground for a principally new scenario of the galaxy for- 
mation in the model of hot Universe. Traditionally, the 
hot Universe model assumes a homogeneous distribu- 
tion of matter on all scales, whereas the appearance of 
observed inhomogeneities is related to the growth of 
small initial density perturbations. However, an analy- 
sis of the cosmological inferences of the theory of ele- 
mentary particles indicates the possible existence of 
strongly inhomogeneous primordial structures in the 
distribution of both the latent mass and baryons. These 
primordial structures represent a new factor in theory of 
galaxy formation. Topological defects such as the cos- 
mological walls and filaments, primordial black holes, 
archioles in the models of axion cold latent mass [19, 
20], and essentially inhomogeneous baryosynthesis 
(leading to the formation of antimatter domains in the 
baryon-asymmetric Universe) [12, 21] offer by no 
means a complete list of possible primary inhomogene- 
ities inferred from the existing elementary particle 
models. 

The proposed approach discloses a number of inter- 
esting aspects in this direction. Indeed, this model pro- 
vides for a possibility of the quantitative analysis of 
correlations in the formation of massive PBHs and the 
primary inhomogeneity of the latent mass and baryons. 
Originally inherent in this mechanism is the inhomoge- 
neous phase distribution which eventually acquires 
(similar to what takes place in the invisible axion cos- 
mology) a dynamical sense of the initial amplitude of 
the coherent oscillations of a scalar field. Irrespective of 
the efficiency of dissipation of the energy of these oscil- 
lations, the regions of closed wall formation must be 
correlated with the regions of maximum energy density 
of the latent mass. If these oscillations are not decaying, 
their energy density may provide for the contemporary 
latent mass density. Inhomogeneity in the initial ampli- 



tude of these oscillations would then imply an inhomo- 
geneity in the initial energy density and, hence, the 
regions of black hole formation will become the 
regions of increased latent mass density. Qualitatively 
similar effect (albeit not as pronounced) takes place in 
the dissipation of coherent oscillations at the expense of 
particle production. An increase in the oscillation 
energy density transforms into a local increase in the 
density of latent mass particles produced in this region. 

In the class of spontaneous baryosynthesis models, 

a change in the phase determines the production of 
excess baryons. Therefore, in addition to an increase in 
the latent mass density, the regions of massive PBH for- 
mation may be characterized by a higher baryon den- 
sity. Inside a closed wall, where the phase is 9 > 71, the 
same mechanism leads to the production of excess anti- 
baryons [12]. However, this antimatter domain will sur- 
vive only provided that its size is sufficiently large [22]. 
Thus, development of the proposed approach may lead 
to a number of interesting scenarios of initial stages in 
the formation of protogalaxies, depending on the selec- 
tion of particular elementary particle models and their 
parameters. This study presents the first step in this 
direction. 
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